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Q ■ Abstract 

We theoretically report that, at a sharp electrostatic step potential in graphene, massless Dirac 
fermions can obtain Goos-Hanchen-like shifts under total internal reflection. Based on these results, 

u ; 

we study the coherent propagation of the quasiparticles along a sharp graphene p-n-p waveguide 
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and derive novel dispersion relations for the guided modes. Consequently, coherent graphene 
devices (e.g. movable mirrors, buffers and memories) induced only by the electric field effect can 
be proposed. 

PACS numbers: 81.05.Uw, 73.63.-b, 73.21.-b, 42.25.-p 



A totally reflected light beam can laterally shifted from the position expected by a ge- 
ometrical trajectory, an effect termed the "Goos-Hanchen" (GH) shift [l|. This effect can 
be qualitatively attributed to evanescent wave penetration into the medium with a smaller 
refractive index. Some theories, such as stationary phase models [2], current-flux models 
[3|], etc., have been developed to formulate this effect. Recently, the exploration of negative 
refraction brought on new interest in this subject since it can lead to negative GH shifts 
[4J, |5|]. Based on negative GH shifts and waveguide theory 6J, the "trapped rainbow" storage 
of light in a tapered left-handed heterostructure (LHH) offers a novel mechanism to bring 
light to a complete standstill efficiently and coherently, which may lead to applications in 
optical data processing, storage and quantum optical memories [7|. 

Graphene, a monolayer of carbon atoms, has attracted great interest since it was success- 
fully prepared [8|, |9|, [lOj. The interaction between electrons and the two-dimensional (2D) 
honeycomb lattice of carbon atoms results in a gapless band structure near Dirac cones with 
a linear (photon-like) electron-hole dispersion, which gives rise to new quasiparticles. Be- 
sides, the hopping of electrons between the two triangular sublattices introduces pseudospin 
1/2 and chirality to the quasiparticles, which can be analogous to massless neutrinos J9l.ll0|. 
Therefore, the quasiparticles can possess the low-energy dynamics effectively described by 
a Dirac equation with Fermi velocity vf ~ 10 6 m/s, and they are also called massless Dirac 
fermions. Furthermore, long phase coherence length 11| and tunable Fermi levels based on 

nn 

electric field effects [8|, [12|] have been experimentally observed in graphene. Based on the 
above properties, recent theoretical research has focused on the coherent transport proper- 
ties of quasiparticles in the electrically gated graphene structures, including Klein tunneling 



13l . [14L focusing of electron flow 15|, confined states and resonant tunneling in quantum 



wells J16l. Il7|. coherent transmission through graphene strips J18l. Il9j]. etc. All these studies 
imply that it might be hard to fully and coherently control the group velocity of the mass- 
less quasiparticles in the graphene structures induced by electric field effects alone. In some 
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221 ] or magnetic fields 



231 ] can be introduced, but for 



cases, graphene quantum dots 
trapping only. 

In this paper, however, we have discovered the possibility to coherently decelerate, trap 
and re-accelerate the massless Dirac fermions along a sharp graphene p-n-p waveguide by 
using a tunable electric field. The waveguide in our scheme can be produced in a bulk 
graphene sheet by electrostatic gates only, which avoids the complex effects of the graphene 
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FIG. 1: (Color online) (a) Total internal reflection and negative GH-like shift at a sharp p-n 
junction, where 9 is the incident angle and d is the lateral displacement of the shift. The gradient 
intensity area represents the evanescent field. On the n side, the pseudospin S of an electron- 
like quasiparticles is parallel to the momentum k\. (b) Low-energy linear spectrum and potential 
diagram. 



strip edges [18 
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Our results are the electronic 



counterpart to the trapped rainbow effect in optics [7|. Before discussing the coherent 
propagation of quasiparticles along a sharp graphene p-n-p waveguide, we first examine the 
total internal reflection (TIR) of quasiparticles at a sharp electrostatic step potential and 
focus on the GH-like shift, which is inspired by the similarities between the reflection and 
refraction of photons and those of graphene quasiparticles at interfaces [10l . [15]. 



In our scheme, we theoretically consider the coherent 



and ignore the disorder as was also done in Refs. 
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231 ] . Moreover, since the electron-electron interaction can not severely distort the shapes of 



constant energy contours of low energy Fermi levels 24j , it can be ignored as well 13| 
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171 . Il8l Il9j . We also assume that, the edge smearing length of the sharp potential is much 



smaller than the incident quasiparticle wavelength, but much larger than the lattice constant 
of graphene, which can prevent intervalley scatterings [l4J. The quasiparticle dynamics is 



therefore governed by the 2D Dirac equation 

[hv F (*-k) + V{x,y)]i/> = Eip, (1) 

where v F ~ 10 6 m/s is the Fermi velocity, cr = (a x , a y ) are the Pauli matrices, k = (k x , k y ) = 
—i'V is the momentum operator, the potential V(x,y) = VqH(x) is y-independent with the 
Heaviside step function H, and E is the eigenenergy. Because a unit cell of the graphene 
honeycomb lattice contains two sublattices A and B, the state ip is expressed by the two- 
component spinor ip = (ip a(x , y) , ip b(x , y)) 7 ', where ipA^x.y) and i(jB{x,y) represent the 
smooth enveloping functions at each sublattice. The conservation of k y leads to ip m {x,y) = 
(x)e lkyV , m = A, B. Thus, Eq. (OQ) gives 
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- ihv F (d<j) B /dx + ky<t> B ) = (E- V H(x))(}>a, (2a) 

-ifiv F {d(j) A /dx - k y (f>A) = (E- V Q H(x))<j) B . (2b) 

The wave functions in the two different regions can be written in terms of incident, reflected 
and evanescent waves with the incident angle 9. Figure [T] shows a p-n junction for clari- 
fication. We define k\ = \E\/hv F , k 2 = \E — V \/hv F , and k y = kism9. In region J, we 
have 

1>i = i/>i + Vi (3) 
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where k x = k\cos8, s = sgn(E), ip\ (ip[) is the incident (reflected) wave function and r is 
the reflection amplitude. In region II, the evanescent wave is 

*n=\(., f ,\ w , )e— + ^y, (4) 



y/2 \is'(a + k y )/k 2 

where the decay constant a = {k 2 — k^) 1 ^ 2 = (k 2 sm 2 9 — fc|) 1//2 , s' = sgn(E — Vq), and t is 

the transmission amplitude. Because a and k x are nonzero real numbers, we can obtain the 

inequalities 

{ -^^<kl<^^0<V <2E, (5) 

n v F nvp 

which is a necessary condition to obtain evanescent waves. When condition (J3]) is satisfied, 
the wavevector k% in region II is smaller than ki in region /, which is analogous to the 
propagation of light from a medium having a larger refractive index into one having a smaller 



refractive index. Hence, the critical angle 9 C for the TIR can be defined by sin# c = A^/fci. 
One can calculate the reflection amplitude r with the boundary conditions of the wave 
functions at x — 0, 



1+r = t, 
ssm9 c (e w - re~ ie ) = is't(6 + sin 9), 



where 5 = y sin 2 9 — sin 2 9 C . From Eq. (jHJ), we can get 

sin 9 C cos 9 + ry 

r = \r\ exp[i(p r \ = ^— a ~ a -, 

sin 9 C cos 9 — i^ 

where 7 = sin0 c sin# — ss'(S + sin 9) and |r| = 1. 
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FIG. 2: (Color online) Negative GH-like shift (solid) at the p-n junction (interband tunneling) 
with s = 1 and s' = — 1. Positive GH-like shift (dashed) at the n-n' junction (intraband tunneling) 



with s 



1. In region I, the wavelength of the quasiparticles is Ai = 2ir/k\ = 100 nm. In 



region II, A2 = 2ir/k2 = 200 nm. The critical angle is 6 C = 30°. 



Some further discussions follow. First, we compare the TIR of the massless Dirac fermions 
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261 ] because both of 



in graphene and that of the Dirac electrons in high energy physics 

them are (pseudo)spin 1/2 fermions. For the TIR of the Dirac electrons, the incident and 

reflected momentum states are well defined due to the conservation of momentum, but the 



spin has two states (e.g. parallel and anti-parallel to the momentum direction), which means 



there are two eigenstates for incidence and reflection, respectively. Different eigenstates have 
different phase shifts, which can cause the splitting. However, for the TIR of the massless 
(neutrino-like) quasiparticles in graphene, the incident and reflected momentum states are 
well defined, and as a consequence, the pseudospin states are well defined due to the chirality. 
This means there is only one eigenstate for incidence and reflection, respectively. Similar to 
Dirac electrons, one can calculate and find the rotation of the pseudospin of the quasiparticles 
can not produce any extra phase shifts and the reflection amplitude r provides the net 
phase shift for the reflected wave. Our analysis indicates that not only does the chirality 
of the massless quasiparticles ensure the absence of back-scattering for normal incidence 
[131 . ll4J . but also no splitting under TIR, such as in Ref. [16(. Second, we also note the 
research on the nonlinear screening in the graphene p-n junction [27[, which can lower the 
junction resistance. Nevertheless, in our scheme [e.g., Fig.QJa)], the shift originates from the 
evanescent wave penetration into the p side having a smaller effective refractive index. The 
nonlinear screening near the interface can not severely affect the existence of the evanescent 
wave on the p side. In principle, this effect can thus be ignored in our theor y 

Based on the above discussions and the stationary phase treatment [J, [25|, [26|] , the lateral 
displacement is 

dtp r 1 d(f r 



d = sm9- 

'r)k x k\ '69 

(8) 



dk x k\ 39 

sin 8 c (ss' cos 2 9 C sin 9 + 5(ss' — sin 9 C )) 



kiS(5 + sin 9) sin 9(1 — ss' sin 9 C ) 
Note that Eq. ([8]) is not valid at the critical angle 9 C and 9 = 90°. The results depicted in 

n n n 

Fig. [2] can be directly analogous to the positive and negative GH shifts in optics [ll, U, l5| . 
Equations ([7]) and ([8]) also tell us that the phase and position of the reflected wave can be 
coherently controlled by adiabatically changing the gate voltages, which can make the step 
potential work as a "movable mirror" . 

The negative GH-like shift of massless Dirac fermions at the sharp p-n junction suggests 
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that, in analogy to the trapped rainbow effect of light in the LHH [7[, the quasiparticles can 
have novel propagation properties in a p-n-p waveguide. For simplicity, we discuss here a 
symmetric p-n-p waveguide [Figs. EJ^a) and[3^b)]. Therefore, Eq. ([2]) is still applicable, but 
the solutions and boundary conditions are different. In region III, we have 

1>iii = a( * V*+^. (9) 

\—is'(o — sm9)/ sm& c J 
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FIG. 3: (Color online) (a) Schematic of a symmetric p-n-p waveguide with a ray in region / and 
evanescent fields in regions II and ///, where the incident angle is sin# = k y /k\ = frvfk y /E. (b) 
Low-energy linear spectrum and potential diagram. 



In region /, we have 



B cos(k x x) + C sin(k x x) 
K is(B sm(k x x + 9) — C cos(k x x + 9)) 



ik y y 



In region II, we have 



1>n = D[ . , 



-a(x—w)+ik y y 



A -- 


= B, 


A(5-sm9) = 


= ss' sin 9 C (B sin 9 — C cos 


D = 


= B cos k x w + Csmk x w, 


D(S + sin9) = 


= ss' sin 9 C (B sin(k x w + 9) 




— C cos(k x w + 9)). 



(10) 



(11) 



Js'(5 + sin 9)/ sin 9 C 
where w is the width of region I. The boundary conditions result in 

(12a) 
(12b) 
(12c) 

(12d) 

A transcendental equation can be derived from the coefficient determinant of Eq. f)12p . which 

is 



sin 9 c (—ss'S cos 9 cos(kiw cos 9) + 

(sin 9 c — ss' sin 2 9) sin(kiiu cos 9)) = 0. 



(13) 



The preceding relationships between k±, &2, E, and Vq give the dispersion relations between 
E and k y for the guided modes. Condition (jSJ) leads to Vq/2 < E < Vb, s — 1, s' = — 1 for the 
p-n-p waveguide. Numerical calculations of the dispersion relations clearly indicate the slow 
positive (forward), zero and even negative (backward) group velocities of the quasiparticles 
depending on the gate voltages [Fig. HJa)] and the widths [Fig. 11(b)] . 

More interestingly, the field-effect-dependent dispersion relations shown in Fig. H(a) 
indicate the possibility to dynamically and coherently decelerate, stop and re-accelerate the 
guided quasiparticles. The detailed process is as follows. First, it is known that, when the 
quasiparticles are guided in the waveguide, the incident angle 9 should be preserved due to 
the parallel interfaces [Fig. [3(a)] . Thus, assuming the quasiparticles are initially prepared 
in a forward propagating state shown at point A in Fig. 11(a), then one can adiabatically 
changing the gate voltages to transfer the forward propagating state to a trapped state at 
point B due to the preservation of the incident angle. After a certain storage time, the 
trapped state can be transferred to either the original state or a backward propagating state 
shown at point C. Consequently, it can be seen that the p-n-p waveguide is not only a 
waveguide, but also a "buffer" or "memory" for the quasiparticles in a coherent way. 

To explain this process, we adopt the current-flux method |3( instead of the stationary 
phase treatment mentioned above. The current density J in the different regions can be 
expressed by 

Jn = VfilpltTlpn) (14) 



with n = I, II, III for the corresponding regions [16J. For the guided modes, we are only 
interested in their propagation along the p-n-p waveguide [the y direction in Fig.[3(a)]. 
Hence, the flux $, which also gives the group velocity v g of the quasiparticles, is given by 

" + 0O 



$ 



/-t-oo 
(Jl,y + Jll,y + JlII,y)dX, (15) 

■oo 



where Ji yV for the guided components gives the positive current density in region /, Ju,y 
and Jin, y for the evanescent components give the negative current densities in regions II 
and III . By adiabatically changing the gate voltages, we can coherently redistribute the 
current densities in the different regions. Consequently, the dynamics of the whole wave 
function are strongly influenced, and the group velocity can be coherently controlled. Our 
calculations show that Eq. (IT5l) gives the same values of group velocities at the points A, B 



and C in Fig.H^a). The current densities in the p-n-p waveguide are schematically illustrated 
in Fig.[5j 

Additionally, the width-dependent dispersion relations shown in Fig. IH(b) also indicate 
the possibility to coherently decelerate and stop the guided quasiparticles by adiabatically 
reducing the width from point D to F while keeping the gate voltages unchanged, which is 
similar to the principle of Ref. 7|. After the storage, one can adiabatically change the gate 
voltages for retrieval. 

In summary, we have demonstrated the GH-like shifts of massless Dirac fermions at an 
electrostatic step potential in graphene and the novel dispersion relations for the guided 
quasiparticle modes in a sharp graphene p-n-p waveguide. Our analysis shows that, based 
on the tunable electric field effects, the electrostatic step potential can function as a movable 
mirror for the incident quasiparticles, and the p-n-p waveguide can coherently buffer, store 
and retrieve the guided quasiparticles. These results may have interesting applications in 
the graphene-based quantum information processing. 
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FIG. 4: (Color online) (a) Dispersion curves for different electrostatic potentials, Vo = 18.84 meV 
(dash-dotted), 19.27 meV (dashed) and 20.72 meV (solid), but the same width of w = 80 nm. For a 
particular incident angle 9 = 13.7° (dotted), one can obtain the different group velocities: 48279.4 
m/s at point A (6.3 ^m -1 , 17.7 meV), zero at B (6.5 /mi -1 , 18.2 meV) and -53399.7 m/s at C (6.8 
[xm. , 18.8 meV) on the corresponding curves, (b) Dispersion curves for different widths, w = 80 
nm (solid) and 60 nm (dashed), but the same electrostatic potential of Vo = 20.72 meV. For a 
particular incident energy E = 20.35 meV (dotted), one can obtain the different group velocities: 
493403.4 m/s at point D (15.6 //m _1 , 20.35 meV) and zero at F (3.7 ^m" 1 , 20.35 meV) on the 
corresponding curves. Note that our parameters only lead to one single dispersion curve for each 
pair of w and Vo, and the solid black squares give the cutoffs for each case. 
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FIG. 5: (Color online) Current density diagrams at the points A, B and C in [Fig. 0(a)], where 
(a) corresponds to A, (b) to B and (c) to C. The arrows schematically show the magnitudes and 
directions of the y components of the current densities in the different regions. 
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